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. Abstract Using chiral supersymmetry, we show that the massless Dirac equation in the Taub-NUT gravita- 
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1. Introduction 

The Dirac equation pijj = in the Kaluza-Klein (KK) monopole field — obtained by imbedding the 
Taub-NUT solution into five-dimensional gravity with trivial time dependence [1] — was studied in the mid 
eighties [2]. It was found in particular that there were no zero modes and therefore the Rubakov effect 
was absent. It was later realized, that the motion of a spin particle in the KK monopole field is exactly 
pi i' soluble, and that it has a Kepler-type dynamical symmetry [3], [4]. Recently [5], a Runge-Lenz vector was 

constructed also in the spinning case. Even more recently van Holten found, by applying the geometric 
technique of Ref. [6], that the Taub-NUT gravitational instanton has an = 4 supersymmetry [7]. In 
this Letter, we point out that the fermion problem is also exactly soluble^ rederive the above symmetries 
' following an entirely different approach, and use them to find the spectrum and the S-matrix. The clue is 

H ; . rk f ■ 

I chiral supersymmetry (SUSY): the Dirac operator, p = \ ^ I is supersymmetric in 4 dimensions; Its 

square, 

(Ll) #^(^*%^,).('^' ^J.i/. 

is in fact the Hamiltonian of a supersymmetric quantum mechanical system, with 'fermion operator' 7^ = 
j and intertwining transformations 



in 
o 



(1.2) U ^ -^T and U'^ = U^ ^T^ ^ 



'Hn V Hn 



Indeed, U"^ HqU = Hi. 
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If ip is an iJo-eigenfunction with eigenvalue > 0, then = T^tp is an ifi-eigenfunction with the 
same eigenvalue. Hi has therefore the same spectrum as Hq up to zero-energy ground states which arise as 

solutions of = 0. If Ho'ip = E'^ij;, then if) ^^^^ = ( ii/f) ' ^° eigenfunctions of the Dirac 

operator can be found from those of p"^. 

Let us assume that one of the partner hamiltonians, say Hq, is simple so that its spectrum and sym- 
metries are known. Then the corresponding quantities for Hi are readily obtained by SUSY. In the same 
spirit, if some operator is conserved for Hq i.e. [Hq, Aq\ = 0, then 

(1.3) Ai = U'^AqU 

is conserved for H\, [Hi,Ai] = 0, and the diagonal operator A = diag(Ai,j4o) is conserved for the Dirac 
operator, [p,A] = 0. (With a shght abuse of notations, Ai will be identified with diag(Ai,0), etc.). The 
bosonic symmetries of p can hence also be obtained from those of Hi and have plainly the same algebraic 
structure as those of Hq. They can be, in principle, calculated using Eq. (1.3). Their exphcit forms are, 
however, rather complicated as well as useless, since all their properties follow from those of Aq. They will 
therefore be omitted. 

The supcrsymmctrics of H are easily identified and contain Qi = ij). Those which commute with Qi 
form the supersymmetry algebra of the Dirac operator. This is precisely what happens for the Taub-NUT 
metric: due to self-duality, one of the partner Hamiltonians reduces to the spin Hamiltonian, and can be 
hence exactly solved. The other partner hamiltonian is complicated, but SUSY dispenses us from solving it. 

2. The Dirac operator in the Taub-NUT background 

The Taub-NUT gravitational instanton is described by the 4-metric 

^ld,,-,r.A.d,^' 

where 



(2.1) v\dr'^ -h r'^{dQ'^ -h sin^ Qd^^)^ -h :^{# - 4to1- df}' 



4/77, _* 

(2.2) V = \^ and ^ • df = cos Qd4>. 

r 

In the usual context, the Taub-NUT parameter, m is positive; m < arises, e.g. in the long-range scattering 
of self-dual SU{2) monopoles [3], [4]. The curved-space gamma matrices, -y^, and the spin connection, T^, 
are expressed as [8] 

(2.3) 7^=^..- ^ h 7^-' ^ 



and 



respectively. The operator 9^ commutes with the Dirac operator p = ^^{d^ + ^r^), since nothing depends 
on tj). Therefore, we can restrict ourselves to a fixed eigensector of d^. Requiring ^ to be periodic with period 
167rm, the eigenvalues of —Amid^, denoted by q and identified with the electric charge, are half-integers. 
For a fixed q (assumed positive for definiteness), the Dirac operator p becomes 



(2.5) 



Tt 
T 



/ '^.i^+J-^' 

^a-nW+^VV 
V Am 
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where tt = — iV + qA. It is easy to check that T and arc each other's adjoint with respect to the Taub- 
NUT volume element Vd^x, as they have to be. Using the (anti)self-duality property VV = —AmB, the 
square of Ip is found to be (1.1) with 



Hi=Ho + 



(2.6) 



1 \2^r2 



7?" + {^fv 

4m 



where L{) = r X li + qv IS the orbital angular momentum [conserved for ffo but not for Hi]. The partner 
hamiltonians Hi and Hq differ hence in a complicated expression, and it is not at all obvious that they have 
the same spectra. SUSY implies however that this is nevertheless true. 

The 'lower' (i.e. 7^ = —1) sector is simple: all spin dependence dropped out. Hq is in fact the the 
Hamiltonian for a spin particle in the KK field [3], [4] (times the unit matrix). Decomposing into radial 
and transverse components, iJo^o = -E^*o leads to the non-relativistic Coulomb-type equation 



(2.7) 



*n = 0. 



r"^ dr^' dr' ' ' ''2m "" 'r ' ^^^m' 
Inserting Lg = ^(-^ + 1)) ^ = 9i 9 + li • • • and setting 

(2.8) *o+="+W>T(^,<^)(j) and = u.{r)Y^{e,<t>) {^-^ 

where the 1^'s are the 'Wu-Yang' monopole harmonics, we find that both radial functions u± become 

(2.9) u{r) = r^e^'^TiiX + L + l,2L + 2, -2ikr), 
where F is the confluent hypergeometric function and 

(2.10) = - {q/Amf and X = -4m[{E'^ /2) - {q/Amf]/k. 

By Eq. (2.9), the wave function vanishes at the origin. No hermiticity problem arises therefore. Bound states 



only arise for E^ > {q/^mf' which requires m < 0. Square integrability requires 
p = 0, 1, . . ., yielding the spectrum 

(2.11) ^" " 8^ ~q^{n- ^Jv? - q^) , n = g + 1, g + 2, . . 

Such a state is labelled by the quantum numbers 

E 

L = q,q + 1,... 
fi = -L,...,L 
s 



-ip+L+iy 



-n 



(2.12) 



energy 

total orbital angular momentum 
third component of L/Q 
eigenvalue of (73 



The n*^ energy level is therefore 2(n^ — g^)-fold degenerate. By supersymmetry, the positive spectrum of the 
superpartner Hi is again (2.11) with the same multiplicities; the eigenstates carry the same labels (2.12). 
Zero modes only arise in the 'upper' sector. Setting 

(gi(g/4Tn)i/, 




-X 
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where x is a two-component Pauli spinor, the equation p^i = reduces to 
(2.14) [a-n-i{^)V]x = 0. 

Applying the flat-space adjoint operator, t? • 7? + i^q/ArnjV, yields 

(2-15) + 2q-^ + (^) (1 + — ) J X = 0, 

which is precisely the zero-mode equation of the D'Hoker-Vinet 'dyon' [9]. For m > 0, there are no zero 
modes. This is consistent with the known results [2]. For m < 0, Eq. (2.15) admits 2q normalizable solutions, 
described, e.g. in Ref. [10]. They have the lowest possible angular momentum, j = q — ^, and are expressed 
as 



^ — — (g— . . . , (q—^). Scattering states arise for > (q/m)^. For a scalar particle, an incoming/outgoing 
state is characterized by its momentum km = k and kout = k' respectively. The scalar 5-matrix reads [3] 

(2.17) So{k'\k) = J2('^L + l)^^--^V^_^jR-\k')R{k)), 

L>g (^ + ?A)! 

where A is as in Eq. (2.10), „ is the rotation matrix and R{k) denotes the rotation which takes the 2;-axis 

to the fc-dircction. In the spinning case [11], for the Hq dynamics, incoming/outgoing states get an extra 
label, namely the third component of the spin, s and s' , respectively. Since s = k ■ a /2, the .S-matrix of Hq 
is simply 



(2.18) S{k', s'\k, s) = So{k'\k) vl!^^{R-\k')R{k)). 

Now if V'o"^""* is any incoming/outgoing wave for the Hq, then ipi = C/(i/)g"^°"*) is one for Hi with the same 
labels. Therefore, the S'-matrix for Hi is again Eq. (2.18), just like for Hq. 



3. Bosonic symmetries 

Hq is conveniently viewed as the Hamiltonian of a spin i particle — but one with gyromagnetic ratio 
5 = 0. Being uncoupled, the spin vector, 

(3.1) So = ia, 
is conserved. SUSY carries this over to Hi: 

(3.2) Si = U^SqU 

is conserved conserved for Hi. (The explicit expression is not illuminating). The diagonal vector S = 
diag(Si, So) is therefore conserved for ^. Note that Si and Sq are separately conserved for H = diag{Hi, Hq). 

The Taub-NUT metric is radially symmetric. Therefore, the orbital angular momentum, Lq, arises nat- 
urally when Ho is viewed as describing two, uncoupled, spinless particles. But the total angular momentum, 

(3.3) J = Lo + So, 
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is also natural, when Hq is viewed as describing a spin^ particle. Hq admits therefore two independent 
conserved 'angular momentum' operators. The operators So and Zq plainly commute, so we get an o(3)®o(3) 
symmetry. Again by supersymmetry, the superpartner, Hi, admits the same symmetries, namely Si in Eq. 
(3.2) and 



(3.4) 



Li = U^LoU. 



Since the total angular momentum, J, is invariant, W JU = J, the partner of Zq = J — So is also written 
as Li = J — Si = Lo + So — Si. The Dirac operator p admits therefore the bosonic o(3) ® o(3) symmetry, 
generated by 



(3.5) 



S = diag(Si, So) and Z = diag(Zi, Lo)- 



These operators were found by van Holten using a quite different method. For example, L is his 'improved 
angular momentum'. 

Hq actually admits further symmetries [3], namely a Runge-Lenz vector, 



(3.6) 



The vector operators Lo and J?o generate an o(3, 1) dynamical symmetry for scattered motions and o(4) for 
bound motions, to which So adds an extra o(3). 

SUSY transports the Runge-Lenz vector to Hi: 



(3.7) 



Ki = U^KqU, 



and then to Ip to yield K = diag(i?i, i?o)- (The explicit form of Ki — presented in Ref. [5] — is very 
complicated). The dynamical symmetry can be used to derive the bound-state spectrum and the S-matrix: 
for fixed energy H = E'^, the Casimirs of the full system satisfy, by construction, the same relations as those 
of Ho, namely 



(3.^ 



K 



1/^2 - (g/4m)2 
K 



■ L 



-4mq 



( 



E^/2-{q/4.mf\-^ 



^E^ - {q/Amy 



confirming that the spin is indeed decoupled. From the representation theory of o(4) we infer that 

. (g/4m)2 - £;2/2 
(3.9) Am^^^==^===L= = n, n = q+l,. 



^J{q/Amf - E'^ 

yielding the spectrum (2.11) once again. The scattering can be described using the spinning extension of 

Zwanzigcr's approach [11]; one finds Eq. (2.18) once again. It is worth stressing that the only effect of spin 
is to double the multiplicity for bound states, and to multiply the S matrix by the spin factor, respectively. 
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4. Supersymmetry [12] 

Consider first the supersymmetries oi H = p'^, restricted to positive-energy states. Firstly, the two odd 
scalar charges 

(4.1) Qi={ rr n ' <52 = -^7 <9i 



TO)' \iT 

are both square-roots of H, 

(4.2) {Qa,Q0} = 25^f}H. 

Now the components of the two 'spins', So and Si, span two independent o(3) ~ su(2)'s, 

(4.3) [SlSl]=zep,,St 
(a, 6 = 0, 1) as well as commute with 7^, 

(4.4) [7M]=0. 

But they also satisfy the anticommutation relation {S„, S^} = (1/2) 6^^6ab', therefore, the following two vector 
odd charges 

(4.5) g„ = 2i[So,Q„] (a = 1,2) 

also commute with H = ij)'^ . The commutation relations are (4.3)-(4.4), augmented with 

[7^ Qo\ = ^iCapQ/}, [7^ Qa] = 2«ec«/3Q^, 

(4.6) [5*, Qi] = li{SijQ^ + Ci^kQi) , [SlQi] = -\i{5ijQo, - ajkQ^) , 

The anticommutation relations read in turn (4.2) and 

{ga,Q/3}--4i/e,^(Si+So), 

Since all odd operators anticommute with 7^, we conclude that H = admits the superalgebra u(2/2) as 
symmetry. It is easy to read off the symmetries of p = Qi from these relations: Qi itself and 

. / ^(^ + i^x^)-(3^)yya\ 

(4.8) 02 = 

\-^{if-iax7r)VV-{^)VVa / 

generate, with S = diag(Si, S2), an A'' = 4 sub-superalgebra: 

[S\Qi]=ieijkQ^^, [S,Qi]=0, 

(4.9) . , . . 
{Qi, Oi} = 2H, {Ql Qi} = 2HSij, {Qi, Q2} = -^HS. 

In Ref. [7], these supercharges are associated to Killing- Yano tensors. 
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5. Discussion 

Some of the results above can be generalized. The multi- Taub-NUT metric is given, for example, again 
by (2.1), but (2.2) replaced rather by 

Requiring (anti) self-duality, VV = —4mB, we still obtain a solution of the vacuum Einstein equations. The 
square of ^ is (1.1) withifo as in Eq. (2.6) (while Hi is a complicated expression similar to the one before). 
The spin drops out again in the 'lower' sector, giving rise to the same supersymmetry algebra as above (while 
the o(4)/o(3, 1) Kcplcr symmetry is in general broken). Similar statements should hold for other self-dual 
metrics like the Eguchi-Hanson or Atiyah-Hitchin metric, etc. 

Unlike in Ref. [5], our conserved electric charge q contains no spin contribution. The explanation may 
well be the following one. Taub-NUT space has the geometry of x R+; the four- dimensional rotation 
group, S0(4), acts on the part by isometrics. The = const, condition breaks S0(4) to S0(3) x S0(2); 
the conserved quantities associated to the S0(3) (respectively to S0(2)) are the total angular momentum, 
J, and the electric charge, q, respectively. The conserved quantities can, however, be 'improved' [7]. For the 
S0(3), this yields the 'orbital' expressions L = J — S. It is likely that the spin part of the 'charge' of Ref. 
[5] can be removed by a similar 'improvement', leaving us with the spin-independent —id^. 

In the 'dyon' case, an additional conserved quantity is provided by the analog of Dirac's operator [10], 

K, commutes with the Dirac operator and is associated to a Penrose-Floyd tensor [13]. It has the remarkable 
property, that its anticommutator with the vector supercharge Q2 reproduces the Runge-Lenz vector. More 
precisely, 

(5.3) \{Q2,lC] = K-q{J + S). 

The results of van Holten indicate, that a similar statement holds in the Taub-NUT case also. We have, 
however, not yet found the curved-space version of the K in Eq. (5.2); the obvious guess ^/VK. neither 
commutes with ]p nor agrees with Eq. (5.3) completely. 

Another intriguing problem is to extend the bosonic o(4, 2) symmetry [3, 4] to the spinning case. 

Acknowledgement . We are indebted to L. Feher for enlightening discussions and for allowing us to use his 
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